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The possibility of the existence of stable tetrahedral deformations at spin zero is investigated 
using the Skyrme-HFBCS approach and the generator coordinate method (GCM). The study is 
limited to nuclei in which the tetrahedral mode has been predicted to be favored on the basis of non 
self-consistent models. Our results indicate that a clear identification of tetrahedral deformations is 
unlikely as they are strongly mixed with the axial octupole mode. However, the excitation energies 
related to the tetrahedral mode are systematically lower than those of the axial octupole mode in 
all the nuclei included in this study. 

PACS numbers: 21.60.Jz, 21.10.Rc, 27.60.+j, 27.50,+c, 27.70.+q 



I. INTRODUCTION 

Exotic shapes of the nuclear density have always at- 
tracted the interest of physicists. In this respect, the 
octupole degree of freedom has played a special role. Ax- 
ial octupole deformations are well established, both ex- 
perimentally and theoretically in several regions of 
the nuclear chart. It has also been shown that non ax- 
ial octupole shapes are competitive with the axial ones 
in specific nuclei 0. However, octupole deformations 
are more subtle than quadrupole ones. Stable static oc- 
tupole deformations correspond usually to shallow min- 
ima as a function of the octupole deformation 0, 0, H| . 
Dynamical studies have shown that octupole correlations 
in the ground state manifest themselves predominantly 
by a spreading of the wave function around the left-right 
symmetric mean- field configuration [||. 

Nuclear tetrahedral deformations have been recently 
investigated in Refs.@, H, Q. Symmetry arguments show 
that tetrahedral deformations induce a 4-fold degener- 
acy of the single-particle spectra. Based on this peculiar 
shell structure, the stability of such configurations has 
been conjectured in specific nuclei. Such a large degen- 
eracy results indeed in large gaps in the single-particle 
spectrum and increases the shell effects for specific val- 
ues of the neutron and proton numbers. Those values 
have been dubbed 'tetrahedral magic numbers' and have 
been determined using a microscopic-macroscopic model 
based on a Woods-Saxon potential @, @, ©, M UjL O • 

However these conjectures almost exclusively originate 
from approaches based on non self-consistent average 
nuclear potentials, with only a limited support of self- 
consistent calculations. It is therefore clearly necessary 
to test their validity in the framework of up-to-date the- 
ories using modern energy density functionals. We ad- 
dress this issue in the present paper in the framework of 
self-consistent mean-field methods using Skyrme interac- 
tions. Since to go beyond a pure mean-field approach has 
been shown to be crucial for octupole deformations, we 
have also studied the stability of tetrahedral shapes us- 
ing the generator coordinate method (GCM). This study 



is focused on nuclei in which the tetrahedral mode has 
been predicted to be favored: 80,98,iio Zrj i 52 -i56 Gd and 
160 Yb. It extends our previous study which was lim- 
ited to 80 ' 98 Zr[H,[i3| (Ref. I) by considering the dynami- 
cal coupling between the axial octupole and the tetrahe- 
dral degrees of freedom. Since tetrahedral shapes are 
generated by the non-axial intrinsic octupole moment 
Q32 oc r 3 (F 3 2 + Yz-2) it is likely that they are in strong 
competition with axial octupole shapes. Most details 
about our method can be found in Ref. I, those on the 
GCM in Ref. [l5| and on its application to 2-dimensional 
octupole calculations in Ref. \u\ . 

Our aim in this study is to determine whether there are 
situations in which a configuration can be unambiguously 
identified as tetrahedral. We will therefore first identify 
which are the possible coexisting structures in the nuclei 
for which tetrahedral deformations have been predicted. 
We will then study whether some of these configurations 
provide clear signatures of tetrahedral shapes. 



II. MEAN-FIELD CALCULATIONS AND 
PARITY PROJECTION 

Octupole deformations of the nuclear density are gen- 
erated by introducing in the HFBCS equations the axial 
Q 30 oc r 3 Y 3Q and the triaxial Q 32 oc r 3 (Y 32 + Y3-2) mo- 
ments as constraining operators. This last one is the op- 
erator generating tetrahedral deformations. The details 
of the HFBCS calculations have been described in Ref. 
I. The pairing interaction strength has been adjusted to 
reproduce 'experimental pairing gaps' in the same way as 
described in Ref. I. In particular the standard prescrip- 
tion based on the odd-even difference of binding energies 
has been used [lsj . The pairing strength for 98 Zr and 
110 Zr has been adjusted to reproduce the pairing gap in 
102 Zr. In Gd isotopes the pairing strength has been ad- 
justed to reproduce pairing gaps in 154 Gd. For studies of 
80 Zr and 160 Yb the pairing strength has been adjusted for 
each of these nuclei individually. These reference nuclei 
have been chosen in order to minimize the influence of 
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deformation changes for the determination of the pairing 
gap. 

The experimental data suggest that all the nuclei in 
which tetrahedral deformations have been predicted the- 
oretically have a well deformed ground state (with the 
possible exception of 98 Zr). A first question to address is 
whether a spherical configuration may coexist at low ex- 
citation energy in any of the studied nuclei. As we shall 
see in the following, the single-particle pattern typical of 
a tetrahedral configuration is very well preserved close 
to sphericity. Therefore a 'tetrahedral magic' nucleus in 
which the spherical configuration is at low energy rep- 
resents a particularly promising case of the existence of 
tetrahedral configurations. Let us first focus on the case 
of 110 Zr for which the variation of the energy as a function 
of the axial quadrupole deformation is shown in Fig. 1. 
Octupole deformations are set to zero for all values of 
the quadrupole moment. The calculations have been per- 
formed for three representative Skyrme parametrizations. 
A deformed prolate ground state at large deformation is 
obtained for all three interactions. It coexists with an 
oblate minimum and, in the case of SLy4 and SkM*, 
with a local minimum for the spherical configuration. It 
is only for the SLy4 parametrization that the spherical 
and the prolate minima are almost degenerate. 

The case of 110 Zr illustrates well the energy depen- 
dence on the quadrupole degree of freedom obtained for 
all studied nuclei with the three Skyrme forces. The 
ground state corresponds in all cases to a deformed pro- 
late minimum (with the only exception of 80 Zr for the 
SLy4 parametrization but this is in clear contradiction 
with the experimental data 16]). The spherical config- 
urations are excited by a few MeV, except in some Zr 
isotopes where depending on the Skyrme parametriza- 
tion they are obtained at excitation energies smaller than 
1 MeV. 




FIG. 1: Variation of the total mean-field energy as a function 
of the quadrupole Q20 moment (in fm 2 ). Triangles, squares 
and circles denote results obtained using the SIII, SkM* and 
SLy4 parameterizations of the Skyrme force, respectively 

A first question that must be answered is how this pic- 
ture is affected by tetrahedral deformations and whether 
the magnitude of the energy gain obtained around the 
spherical minimum is large enough to bring its energy 



close to that of the deformed ground state. 

The variation of the mean-field energy with octupole 
deformations when the quadrupole moment is con- 
strained to be zero is shown in Fig. 2 for the six nuclei 
that we have selected. These results have been obtained 
with the SLy4 Skyrme parametrization and their depen- 
dence on the choice of the parametrization is illustrated 
by results obtained for 110 Zr with SIII. Note that, con- 
trary to Fig. 1, the HFBCS results are very similar for 
both interactions. We have also checked for several nuclei 
and for SkM* that the behaviour of the HFBCS energy 
as a function of the octupole degrees of freedom is qual- 
itatively similar in all cases. 

The octupole deformations are parameterized by di- 
mensionless parameters proportional to the octupole mo- 
ments (see Ref. I). There is no direct connection be- 
tween the values of /3sq and P32 and one must be cau- 
tious when comparing axial and tetrahedral deformation 
energy curves. The only nucleus for which a sizeable 
tetrahedral minimum is obtained is 156 Gd with an en- 
ergy gain of 500 keV with respect to the spherical config- 
uration. The energy is rather flat as a function of both 
/?30 and P32 for the other nuclei, although the meaning of 
this shallowness can only be made precise by dynamical 
calculations. The energy curve obtained for 110 Zr with 
SIII is only marginally different from that of SLy4, with a 
tetrahedral minimum 100 keV below the spherical point. 

These mean-field calculations have been performed us- 
ing a single constraint on either the Q30 or Q32 mo- 
ments. Since the unconstrained degrees of freedom are 
completely relaxed (except for the quadrupole moment 
and either Q30 or Q32 which were set to zero) by the 
variational procedure, there is no guarantee that the un- 
constrained moments remain equal to zero (except those 
that are forbidden due to the imposed symmetry condi- 
tions - see Ref. I) for all values of the constraint. Never- 
theless, the behavior typical of the tetrahedral symmetry 
is largely preserved up to a deformation of ^32 ~ 0.2—0.3. 
Up to these values, the single-particle states exhibits the 
4- fold degeneracies characteristic for the point group . 
In general, the single particle energies as a function of Q32 
exhibit more bunching as compared to the Q30 direction. 
However, this does not translate into a lower energy for 
the tetrahedral configuration. 

As a typical example, the variation of the single- 
particle energies as a function of octupole deformations 
is shown in Fig. 3 for 110 Zr. One can see that the single 
particle states are still almost degenerate at ,832 equal 
to 0.15 but not at P30 equal to 0.15, which is in both 
cases the deformation beyond which the mean-field en- 
ergy starts to increase. Note also that both the tetra- 
hedral and the spherical configurations have a similar 
single-particle structure since there are no level crossings 
between these configurations. The same is true for the 
axial octupole configuration. Moreover the HFBCS cal- 
culations indicate almost no barrier between tetrahedral 
and axial octupole minima. These facts indicate that 
all the three configurations may be strongly mixed when 
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the octupole collective dynamics is taken into account 
(see next section). 

Since tetrahedral deformations break parity, projection 
on parity gives rise to an energy gain for the positive 
parity as soon as the octupole moments have a non zero 
value and it generates a distinct energy curve for the 
negative parity. We have restored both particle number 
and parity by projecting the mean-field wave functions. 
The projected potential energy is defined as: 



E(N,Z,f3 3 ^)± = 



W^)\HP { ±,N,z)\<f>((i^)) 
m^)\P(±,N, Z) m^)) ' 



(1) 



where \<t>{(i 3tl )) are HFBCS wave functions generated 
with the constraint (<M/33m)|Q3mI < M#v)) = C^fo^, where 



Co = j- A 



,C 2 = Co/V2 with A = 



N + Z and 

r = 1.2/to (see Ref. I). The operator P(±,jv,z) is the 
product of operators projecting on tt = ±1 parity and on 
N neutrons and Z protons. The parity-projected ener- 
gies are shown in Fig. 4, except for 80 Zr and 98 Zr which 
were already discussed in Ref. I. 

The situation is similar to that discussed already in 
Ref. [HI, [HI . The positive parity curves exhibit a small 
minimum for nonzero values of /?3o and P32 ■ The energy 
gain due to the parity restoration is of the order of 1 MeV 
for both axial octupole and tetrahedral deformations (see 
Fig. 4). In fact for all nuclei considered (with the excep- 
tion of 160 Yb ) the axial octupole minimum has a slightly 
lower energy than the tetrahedral one (see table I). 

The dependence of energy on octupole deformations is 
not significantly modified by the projection on positive 
parity. An interesting property of the particle projection 
is that it makes the results rather weakly dependent on 
the pairing strength, in contrast to pure mean-field calcu- 
lations. We have checked that the differences between the 
energy of the octupole minima and that of the spherical 
configuration are not significantly modified by a varia- 
tion of the pairing strength within the interval between 
half and twice the physical value. This holds for both 
the positive and negative parity states. The main effect 
of an increase of the pairing strength is an increase of 
the energy difference between the positive and negative 
parity minima. 



III. TWO-DIMENSIONAL GCM 

To be unambiguously identified experimentally, tetra- 
hedral deformations should have a clear signature, which 
allows to distinguish them from axial octupole deforma- 
tions. The GCM enables to study the coupling between 
both octupole modes and to see whether tetrahedral 
shapes can be separated from axial octupole shapes. We 
have therefore performed two-dimensional GCM calcula- 
tions in which the axial and non axial octupole shapes 
are coupled. This coupling was not considered in Ref. 
I, where separate dynamical calculations were performed 



along the collective paths determined by non zero Q30 
and Q32 values, respectively. The method we apply is 
similar to that introduced by Skalski et al. [13] . 

A collective wave function is constructed by mixing 
the mean-field states corresponding to different values of 
the octupole moments, after their projection on particle 
numbers: 

|*>= f f(Paa,fiu)P(N,z) W/%o, /%2)>4&odfta (2) 



The coefficients /(/?3o, P32) are determined by minimizing 
the total energy of the collective wave function \^). 

Our collective space forms a plane specified by Q30 
and Q32, or equivalently by /33o,/?32- This requires to 
consider HFBCS states inside a rectangle specified by 
"corners": (iftoma^ifemaj;)- However the full prob- 
lem can be decomposed in four subspaces by introducing 
combinations of states in the four quadrants. Starting 
from \<fii) = \4>(f33o, P32)), one constructs the four states: 

102) = |0(-/?3O,-& 2 ))=P|0l), 

103) = |0(-/33O,+/?32))=P*yP|0l), 
|0 4 ) = |0(+/3 3O ,-A52))=JW|01), 

where P is the parity operator and P xy is the reflection 
operation in which x and y coordinates are exchanged. 
Thus, one needs to generate only the HFBCS basis in 1/4 
of the rectangle and extends it to the full square thanks 
to these relations. Another interest of this decomposition 
is that P and P xy commute with the Hamiltonian. This 
means that they can be used to label GCM eigenstates. 
Both P and P xy are projectors, so the quantum numbers 
associated to each operator take the values ±1. From the 
wave functions \4>i), i = 1,2,3,4, one can define a new 
basis in which both P and P xy are diagonal, namely: 



*++) 


= \m 


+ |0 2 ) 


+ 103) 


+ 104)), 


$__) 


= \m 


"102) 


+ 103) 


-104)), 


*-+) 


= \m 


-102) 


-103) 


+ 104)), 


*+-) 


= \m 


+ 102) 


-103) 


-104)), 



where the first index of \$kl) denotes the eigenvalue with 
respect to parity and the second index with respect to 

x — y reflection. One can easily check that |$_| ) is 

identically zero in the absence of either axial or tetra- 
hedral deformations, while |<E> |_) and |$__) are zero 

when /?3o or P32 are zero, respectively. For this rea- 
son, we have dubbed the excited states corresponding to 
k — I — — 1 tetrahedral excitations, those corresponding 
to k = — 1, 1 = +1 axial excitations and those correspond- 
ing to k = 1, 1 — — 1 mixed octupole excitations. 

In this basis \$>ki), the Hamiltonian does not couple 
states corresponding to different values of k and I and 
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the GCM equation decomposes into four equations for 
each set (k,l). The resulting GCM wave functions are 
expressed by: 

|*W>= j f(^^2)P(N,Z)W^,^2)kl)d(i^d(i Z2 (3) 

We have restricted this study by imposing the 
quadrupole moment to be fixed. A full calculation would 
require to consider the octupole and quadrupole modes 
simultaneously, which would be a huge computational 
task, well beyond the scope of the present study. Our 
aim is indeed only to determine the most favourable sce- 
nario of coupling between axial octupole mode and the 
non axial Q32 mode generating the tetrahedral deforma- 
tion. 

We have first performed a calculation in the vicinity 
of the deformed ground state of two nuclei, 110 Zr and 
154 Gd. The results are shown in Table II. 

The correlation energy due to the octupole modes is 
defined by: 

Ecorr = E(N, Z,S P h)-E ++ , 

where E(N, Z 1 sph) is the energy of the particle number 
projected spherical configuration and is the lowest 
positive-parity energy obtained in the GCM. The value 
of this correlation energy is rather small for both nuclei, 
indicating a weak effect of octupole correlations in the 
ground state. 

The lowest octupole excitation corresponds in both 
cases to the axial octupole mode. The non axial octupole 
excitation is only slightly larger in energy for 110 Zr but 
both modes being above 4 MeV of excitation are very un- 
likely to survive to the coupling to any other modes. The 
situation is slightly more favorable in 154 Gd, although in 
this case, the non axial excitation is nearly 1 MeV above 
the axial octupole one. Note that in the vicinity of a 
deformed ground state, one cannot identify a non ax- 
ial Q32 mode with tetrahedral deformations, since the 
tetrahedral symmetry is broken by large quadrupole de- 
formations. 

In view of the very unfavorable conditions obtained 
when the quadrupole moment is large, we have continued 
this study by looking to the octupole properties around 
the spherical configuration. The GCM results are sum- 
marized in Table III. We have performed calculations 
with two sets of mean-field wave functions correspond- 
ing to 16 and 25 positive octupole deformations respec- 
tively, to check the accuracy of the results. The difference 
between both sets of results shows that the accuracy ob- 
tained with a 25 wave-function basis set is better than 
100 keV. Note that the excellent agreement obtained by 
using two different basis sets is also a test that our results 
are not affected by the pathology that can appear when 
working an energy density functional (l9j. 

The largest gain is obtained for 110 Zr. It is of a sim- 
ilar order of magnitude than the energy gain due to 
quadrupole correlations in deformed nuclei [20j. A full 



study of the energy gain due to the coupling between dif- 
ferent modes remain to be done but the results of Ref [2(| 
seem to indicate that these energy gains quickly saturate 
in models based on self-consistent mean-field wave func- 
tions. 

Dynamical deformations associated with the lowest 
GCM solutions corresponding to quantum numbers k and 
I are defined by: 

&/. = 2/Wl(03O,0S2)> (4) 

for /1 = 0, 2, where gu is the collective wave function for 
parity k and for an eigenvalue / associated with the op- 
erator P xy (see Ref. [l5[ for the relation between the col- 
lective wave function g and the GCM function /). For all 
nuclei that we have studied the dynamical deformations 
Pso and P32 of the lowest positive parity GCM solutions 
are smaller than 0.1. The ground state collective wave 
function is rather isotropic as a function of Q30 and Q32 . 
It shows a similar spreading as a function of axial and 
tetrahedral octupole deformations. 

The first negative parity state has an excitation en- 
ergy comprised between 1.0 and 2.3 MeV, the tetrahedral 
mode being systematically the lowest one. The largest 
differences — between both octupole modes oc- 
cur for 110 Zr where it is around 0.8 MeV and for 156 Gd 
where it is around 0.5 MeV. 

The ratio between the B(E3) values obtained for both 
modes are given in Table IV. Better than the absolute 
values of these quantities which are not well defined in 
an angular momentum unprojected model, these ratios 
are good indicators whether these states have a specific 
signature in their deexcitation spectrum. In the first col- 
umn are given the ratios corresponding to the transitions 
from the tetrahedral and the axial excited states to the 
ground state. The second column corresponds to the ra- 
tios of the transitions between the mixed octupole states 
to the tetrahedral and the axial excitations. This ratio 
oscillates in all cases around 1 which suggests that the 
spectrum of GCM excitations resembles to a large extent 
a harmonic spectrum. The only noticeably deviation oc- 
curs in the case of 110 Zr where the transition from the 
tetrahedral state to the lowest GCM state is decreased by 
about 30% as compared to the deexcitation of the axial 
octupole vibration. 



IV. CONCLUSIONS 

We have investigated the possible existence of stable 
tetrahedral configurations in nuclei in which they were 
predicted on the basis of non self-consistent models. Our 
calculations have been based on several parametrization 
of the Skyrme interaction, with only marginal differences 
between the results. The coupling between the axial and 
tetrahedral octupole modes has been studied with the 
GCM, in the absence of quadrupole deformations. Our 
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results do not support the prediction that tetrahedral 
deformations should have a definite signature: 

• The susceptibility of the spherical configuration to- 
wards tetrahedral deformations is rather weak and 
pairing effects wash out the shell effects. Moreover 
the tetrahedral minimum is accompanied by an ax- 
ial octupole minimum of similar depth. 

• The correlation energy associated with shape fluc- 
tuations and parity restoration lowers substantially 
the mean-field energy. However the dynamic oc- 
tupole deformations in the ground state state is 
rather small. Moreover axial and non axial oc- 
tupole deformations are strongly coupled. 

• The excitation energies of states associated with 
tetrahedral shapes are systematically lower than 
those corresponding to the axial octupole mode. 
However the B(E3) ratios do not distinguish be- 
tween these modes. 

The prospects for the experimental detection of the 
tetrahedral configurations at spin zero are thus rather 
poor, ft seems that the increased shell effects due to the 
tetrahedral mode do not provide a sufficient condition 
for the existence of a stable tetrahedral deformation. At 
spin zero, stable tetrahedral configurations seem unlikely. 



Their trace may be manifested in nuclear vibrations in 
negative parity bands but the B(E3) values indicate that 
there is no way to distinguish the tetrahedral modes from 
the axial octupole modes by looking to the decay proba- 
bilities, ft should be noted that our study does not rule 
out the possibility of the existence of rotating tetrahedral 
configurations. Several of the predicted tetrahedral nu- 
clei are however strongly deformed in their ground state 
and a mixing of octupole and quadrupole deformations 
would make still more problematic the extraction of a 
tetrahedral signature. 
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FIG. 2: Total mean-field energy obtained with the HF 
(squares) and HFBCS (circles) methods as a function of the 
dimensionless fim (filled symbols) and fi^2 (open symbols) de- 
formation parameters. 



SLy4 SHI 

Sly4 (3 3 o P32 (330 

110 Zr -1.25 -0.98 -1.12 -0.96 
52 Gd -1.09 -0.93 - 
54 Gd -1.13 -0.92 - 
56 Gd -1.16 -1.14 - 
60 Yb -1.19 -1.20 - 



TABLE I: Energies (in MeV) of the positive parity configu- 
rations with respect to the spherical configurations for axial 
(/?3o) and tetrahedral (^32) configurations. 
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FIG. 3: Single-particle energies as a function of octupole 
deformation P30 and P32 for 110 Zr calculated for SLy4 force. 
The positive and negative parity levels are denoted by solid 
and dashed lines, respectively. 



SLy4 E exc (MeV) E corr (MeV) 71" Ttxy /?30 /?32 



J Zr 1.222 +1 +1 0.08 0.06 

4.535 - -1 -1 0.07 0.17 

4.282 - -1 +1 0.15 0.04 

7.423 +1 -1 0.19 0.24 

Gd 2.228 +1 +1 0.05 0.06 

2.892 - -1 -1 0.06 0.10 

1.998 - -1 +1 0.12 0.04 

5.005 - +1 -1 0.12 0.10 



TABLE II: Excitation energies, correlation energies and dy- 
namical deformations of the lowest four states obtained in 
2-dim GCM. tt and ir xy denote the parity and P xy quantum 
numbers, respectively. 
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FIG. 4: Parity and particle number projected energies as a 
function of the octupole deformations /33o and ,832 . In the two 
top subfigures (a and b) corresponding to 110 Zr, the results 
obtained with the SLy4 (circles) and SIII (squares) are com- 
pared. In the other subfigures, only calculations performed 
with SLy4 are shown. The circles and squares denote then 
the energies as a function of (3zo and respectively. In 
all subfigures filled and open symbols refer to positive and 
negative parity, respectively. 



SHI E exc (MeV) E corr (MeV) A 16/25 n n xy /3 3 o 2E p qp (MeV) 2E% P (MeV) 

80 Zr 2.116 0.15 +1 +1 0.07 0.06 2.786 3.518 

2.832 - 0.06 -1 -1 0.05 0.18 

2.854 - 0.01 -1 +1 0.14 0.00 

8.275 - 0.16 +1 -1 0.13 0.17 

98 Zr 0.0 1.184 0.02 +1 +1 0.07 0.04 2.114 2.66 

2.128 - 0.12 -1 -1 0.04 0.25 

1.732 - 0.09 -1 +1 0.19 0.02 

6.628 - 0.12 +1 -1 0.17 0.23 

SLy4 E exc (MeV) E corr (MeV) A 16/25 n n xy /3 3 o hi 2E p qp (MeV) 2E% P (MeV) 

98 Zr 2.660 0.07 +1 +1 0.10 0.08 1.96 1.78 

2.393 - 0.06 -1 -1 0.08 0.21 

2.639 - 0.05 -1 +1 0.18 0.06 

6.127 0.02 +1 -1 0.17 0.17 

3.303 0.01 +1 +1 0.09 0.10 1.612 2.72 

0.01 -1 -1 0.06 0.22 

0.01 -1 +1 0.17 0.06 

0.01 +1 -1 0.16 0.20 

2.791 0.00 +1 +1 0.05 0.06 2.884 2.78 

0.00 -1 -1 0.04 0.13 

0.00 -1 +1 0.11 0.05 

- 0.01 +1 -1 0.12 0.12 

3.054 0.00 +1 +1 0.06 0.07 2.566 3.0 

0.01 -1 -1 0.05 0.14 

0.01 -1 +1 0.12 0.05 

0.00 +1 -1 0.11 0.02 

3.085 0.10 +1 +1 0.06 0.08 2.008 2.742 

0.00 -1 -1 0.05 0.15 

0.06 -1 +1 0.12 0.05 

- 0.02 +1 -1 0.11 0.13 

3.085 0.00 +1 +1 0.06 0.06 3.438 3.06 

0.00 -1 -1 0.05 0.13 

0.02 -1 +1 0.11 0.05 

0.02 +1 -1 0.10 0.12 



110 Zr 
1.764 
2.188 
4.936 



Gd 
2.018 
2.233 
4.922 



Gd 
1.507 
1.857 
4.134 



Gd 
1.072 
1.507 
3.329 



i60 Yb 

1.629 
1.858 
3.893 



TABLE III: Excitation energies, correlation energies and dy- 
namical deformations of the lowest four states obtained in 2- 
dim GCM. Ai 6 /25 denotes the difference in energies between 
calculations performed with 16 and 25 mean-field states, n 
and n xy denote the parity and P xy quantum numbers, respec- 
tively. In the last two columns the two-quasiparticle excita- 
tion energies (neutron and proton) of the spherical configura- 
tion are given. 



SLy4 a b 

9S Zr 0.78 1.05 
110 Zr 0.67 0.72 
152 Gd 0.83 0.90 
154 Gd 0.97 0.91 
156 Gd 1.15 0.89 
160 Yb 1.04 1.12 



TABLE IV: Ratios of the B(E3) values obtained for the tran- 
sitions between the four lowest GCM states. In the column 
denoted by a, the ratio is taken between |3> — ) — > | ( 3>++), 
and |$-+) — > | <&++}, and in b for: |$+-> —> |$— ) and 

!*+->-> !*-+>. 
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